Saddles in the energy landscape probed by supercooled liquids 
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We numerically investigate the supercooled dynamics of two simple model liquids exploiting the 
partition of the multi-dimension configuration space in basins of attraction of the stationary points 
(inherent saddles) of the potential energy surface. We find that the inherent saddles order and 
potential energy are well defined functions of the temperature T. Moreover, decreasing T, the 
saddle order vanishes at the same temperature (T^fCT) where the inverse diffusivity appears to 
diverge as a power law. This allows a topological interpretation of T^^^: it marks the transition 
from a dynamics between basins of saddles (T > T^ct) a dynamics between basins of minima 

(r < Tj^f^j,). 
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PACS Numbers : 61.43. Fs, 64.70.Pf, 61.20.Ja 

The study of the properties of the free energy land- 
scape or/and of the potential energy surface (PES) in 
disordered systems is a topic of current research ^ . The 
universality in the dynamic of systems as different as dis- 
ordered spin-glasses and structural-glasses |^ supports 
the possibility that some universal features at the land- 
scape level control the slow-dynamics in these systems 
1^. Along this line, several recent works have attempted 
to connect both dynamics and thermodynamics of glass- 
forming liquids to landscape properties. 

Numerical investigations of the sampled configuration 
space have been performed for several models of liq- 
uids An important outcome of these studies is 
the demonstration that on cooling the system populates 
basins of the PES associated to local minima (the so- 
called inherent structure IS ||l^) of deeper and deeper 
depth Ql . The number of distinct basins with the same 
depth in bulk systems has also been evaluated . These 
information have been incorporated in a detailed de- 
scription of the thermodynamics of supercooled liquids 

iiiJil- 

In the landscape framework, the dynamics of the sys- 
tem in configuration space is conceptually decomposed in 
a " fast" oscillatory motion (dynamics within a basin) and 
a slow diffusive motion (dynamics among different PES 
basins). Quantitative calculation of the diffusion coef- 
ficient D based on landscape properties have been for- 
mulated within the Instantaneous Normal Mode (INM) 
theory jl2|. The INM approach focus on the properties 
of the local curvature of the PES sampled by the liquid, 
calculated diagonalizing the Hessian (Ti.) matrix of the 
potential energy V. Analysis of the resulting eigenvalues 
and eigenvectors p3[ | allows to evaluate the number of 
independent directions in configuration space associated 
to basin changes, i.e. to diffusion. For the cases where 
such analysis has been performed, strong evidence has 
been presented for existence of proportionality between 
D and the number of diffusive directions fl^-p^. It has 
also been shown that the number of diffusive directions 



decreases with T and appears to vanish at the so-called 
Mode-Couphng ||l^ transition temperature T^^^ Jl5| , p7t , 
i.e. at the T where an apparent divergence of the inverse 
diffusivity is observed. 

A major difficulty in quantitatively pursuing the idea 
of a connection between diffusivity and topology of con- 
figuration space is that no transparent mapping has been 
yet proposed to associate equilibrium configurations to 
the "closest" configuration on the border between differ- 
ent basins. 

In this Letter we propose such a mapping and present 
an analysis of the properties of the closest stationary 
points (saddles) of the potential energy. The proposed 
mapping, which can be considered as an extension of the 
Stillinger- Weber mapping jl^Jl^, partitions the configu- 
ration space R'^^ of a 3-dimensional TV particle system 
in basins of attraction of the saddles, characterizing the 
saddle with their order and their potential energy. The 
dynamics of the system is then described as dynamics be- 
tween different saddles' basins. We apply the proposed 
mapping to two different models of simple liquids and 
find that: (i) the order of the sampled saddle is a well 
defined function of T; (ii) on cooling the liquid populates 
basins associated to saddles of lower and lower order; (iii) 
the location in potential energy of the saddles is much 
smaller than the system potential energy, providing evi- 
dence that the diffusion process is entropy driven, even 
below Tj^,^^; (iv) at r^j^^, the saddle order appears to 
vanish, indicating that, at this T, the system populates 
the basins of potential energy minima, confirming that 
Tj,^py marks the crossover between two different dynam- 
ical processes. 

To partition R'^^ in basins of saddles we search for the 
basins of attraction of an auxiliary potential function, 
namely W = ^\VV\^ (for a similar approach see also 
^l|). The function W is never negative and it is zero at 
all saddle points, i.e. at all points where all forces are 
zero (stationary point configurations). The saddle points 
are classified according to their order Ug (the number of 
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negative eigenvalues of W ) and their potential energy 
Cs. Saddles of order zero coincides with local minima 
of the PES (i.e. with IS). The complete description of 
the energy landscape would require the calculation of the 
densities of states for each n. However, by investigating 
two model systems we find that similarly to the case of 
minima the saddle's energy and order are well de- 
fined function of T. Hence, all relevant information are 
contained in the functions es{T) and ngiT). 

We investigate numerically two simple model liquids: 
(a) the monatomic Modified-Lennard-Jones (MLJ) [ p2| , 
and, (&) the standard Lennard- Jones 80/20 binary mix- 
ture (BMLJ) Both models are able to support 
strong supercooling without the occurrence of crystal- 
lization. Standard LJ units are used hereafter. Equi- 
librium configurations are prepared by standard micro- 
canonical molecular dynamics simulations at constant 
density (p = 1 for MLJ and p = 1.2 for BMLJ) and at Ts 
ranging from the normal liquid phase (T w 1.6), down 
to r„cT- {Tmct = 0-475 for MLJ and 0.435 for BMLJ). 
The systems are composed of =256 (MLJ) and iV^lOOO 
(BMLJ) particles enclosed in a cubic box with periodic 
boundary condition. Truncated {Rc = 2.6 and 2.5 re- 
spectively) and shifted LJ potentials are used. We an- 
alyze 20 independent equilibrium configuration for each 
T. For each configuration we calculate the associated IS 
and inherent saddle implementing a steepest descent al- 
gorithm which moves in the direction of —'S/V = F and 
—VW = Ti ■ F respectively (the arrows indicates 3iV- 
d vectors). Finally, the H of the starting equilibrium 
configuration (to evaluate the INM) and of the inherent 
saddle (to evaluate Ug) is calculated and diagonalised. 

Fig. I^a shows the average order of the inherent sad- 
dle (i.e. the number of negative eigenvalues of Ti.) as a 
function of the equilibrium T. The average order ns{T) is 
a well defined function of T [|4j , indicating that the tra- 
jectory of the system in the configuration space samples 
statistically the subspace set up by basins of saddles of 
a given order ns{T). Fig. 0a shows also that, decreasing 
T, ns(T) vanishes at T^jct- 

This finding gives support to the following scenario: 
T^MCT is the T above which the system explores basins 
of saddles of order > and below which the system 
is mostly confined in a local minimum (n^ = 0). The 
existence of a quantity, ns{T), that vanishes at T^jp^ is 
remarkable. Indeed, this makes ns{T) a good candidate 
for the description of supercooled dynamics, for instance 
the computation of Us (T) is an alternative way of deter- 
mining T^jcT- 

For comparison with the previous INM studies. Fig. 0b 
also shows the number of directions characterized by neg- 
ative eigenvalues function of T. As discussed in 
Ref. JlSf , a non zero value of n,; is found at T^^j^^, when 
the system is trapped in basins of minima, a clear sig- 
nature of the presence of non-diffusive unstable modes. 
Hence, the introduction of the inherent saddle concept 



offers a way to overcome the difficulties associated to the 
presence of non-diffusive modes in the standard INM ap- 
proach. The order of the saddle Ug appears to be a well 
defined indicator of the number of diffusive directions. 

Next we discuss the location in energy of the inherent 
saddles as a function of the equilibrium T. Fig. || re- 
ports the average instantaneous potential energy e^, the 
average potential energy of the saddle and the aver- 
age potential energy of the IS Bq, as a function of T for 
the MLJ potential model. Similar results holds for the 
BMLJ case. The quantity eo(T) shows a rapid increase 
between Tj^^,, and T w 0.8 reaching a constant value for 
higher T (see inset of Fig. ^ where is reported in an 
expanded scale). We notice that the overall variation of 
Eo is very small on the scale of the variation of and 
Ci. The quantity es(T) shows an intermediate behavior 
between and Co- In agreement with the observation 
that around Tj,,,^^ the system explores basin of attrac- 
tion of saddles of order zero, we find that at T ^^a/ct 
the inherent saddle energy curve merges on the IS curve. 
The data reported in Fig.s |^ and ^ allow to conclude 
in an unambiguous way that T^^^^ marks the cross-over 
between two dynamic regimes: at T < T^^^^ the sys- 
tem spend most of the time trapped in a local minima, 
while at T > T^^^^^ the system explores basins pertaining 
to saddle points of increasingly higher -but well defined- 
energy and order. This conclusion is consistent with the 
INM finding that at T^^^t fraction of diffusive direc- 
tions explored by the liquid goes to zero |jl^,|l5|,0. It is 
also consistent with the interpretation of the transition 
between different dynamical regimes at Tjv/cr proposed 
in Ref. and based on the analysis of the decay of the 
density-density correlation functions evaluated along an 
inherent structure trajectory. 

An important observation stems from the data re- 
ported in Fig.s |l| and |[ the saddles energies are located 
well below the instantaneous values. This finding clearly 
shows that the system trajectory is never close to a saddle 
point and dynamics should not be described as saddle-to- 
saddle motion, but, more correctly as dynamics between 
basins of attractions of the corresponding saddles. In 
this respect, one should not discuss the dynamics of the 
system for T > T^^^^ as activated dynamics It is 

worth to note that the value at T^^^t larger than 
es(T) for a wide T range. In other words, even at Tmctj 
the instantaneous potential energy is much larger than 
the energy value at which saddles are located. Neverthe- 
less at 'rj,^cT tlis system spends a very large fraction of 
time in a local minimum and only rarely performs jumps 
between minima. Hence, the diffusion events at low T 
are not limited by the presence of energy barrier that 
must be overcome by thermal activated processes; they 
are rather controlled by the limited number of directions 
leading from a basin of a minimum to another basin at 
(almost) constant potential energy. 

Fig. 0a shows the relation between and the eleva- 
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tion (in potential energy) respect to the corresponding 
local minima {eg — Co). We discover a remarkable lin- 
ear relationship between these two quantities. This indi- 
cates that -given a minimum- the energy landscape above 
it is organized in families of equally spaced energy sad- 
dle points (to rise one step in the saddle order, the re- 
quested energy is 3NA{es — eo)/Ans ~ 11 for both sys- 
tems. Moreover, this simply linear dependence suggests 
that the aspects of the landscape above a local minimum 
is independent from the energy of the minimum itself. 
We also find a linear relation between the mean square 
distance between minima that have been reached by 
quenching the inherent saddles's of order rig and the or- 
der of the corresponding saddle as shown in Fig. ||b 
This linear relation suggests that the descent path from 
a saddle of order Hs towards the surrounding minima 
can be interpreted as a sequence of independent random 
steps, each of them decreasing the order of the saddle 
by one and increasing the (squared) distance between 
the associated local minima by a fixed amount (random 
walk among saddle points). It is worth to note that these 
properties of the energy landscape (i.e. the linear depen- 
dence between — Co, d^, and Us) are ingredients of a 
model for the landscape introduced by Madan and Keyes 
p7t and recently revisited |28|. This similarity deserves 
further investigation. 

In conclusion, by the numerical investigation of two dif- 
ferent systems, we have pointed out the relevance of the 
concept of inherent saddle in describing the dynamics of 
supercooled liquids and we have highlighted same impor- 
tant characteristics of the energy landscape. In particular 
we have: i) shown that the inherent saddles' properties, 
energies and orders, are well defined function of T; ii) 
demonstrated the validity of the conjecture that T^j^^ 
marks the transition between a dynamic among minima 
(T < Tj,jpy) and a dynamics where the systems spend 
time nearby high order stationary points; Hi) found that 
ns{T) is a quantity that can be efficiently used to mea- 
sure numerically 7^^^^; iv) provided evidence that the 
diffusion processes are entropy driven; and v) shown that 
the aspect of the energy landscape "seen" by a given lo- 
cal minimum is highly regular (as demonstrated by the 
linear dependence of Cg — Co and of c?^ from Us) and in- 
dependent from the specific minimum. 
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FIG. 1. Temperature dependence of the fraction of the neg- 
ative eigenvalues of the Hessian calculated at the inherent 
saddle configurations Us/SN (a), and at the instantaneous 
configurations, Ui/SN (b). Open symbols refer to MLJ and 
full ones to BMLJ. The dashed line in a) are the best fit to 
the data with the function A(T - T^y (MLJ: T„=0.48±0.01, 
7=0.78±0.02. BMLJ: T^=0.435, 7=0.94±0.1.). 

FIG. 2. a) Temperature dependence of the instantaneous 
energy [d, o), the inherent saddle energy [es, o), and the in- 
herent minima energy [co, □) for the MLJ potential model. 
The vertical dashed line indicates T^^^^ as derived from 
the apparent divergence of the inverse diffusivity. The in- 
set shows, in an enlarged energy scale, the T dependence of 
Co- A very similar behavior is found for the BMLJ system 
(not reported here for clarity reasons). 

FIG. 3. a) Saddle's elevation energy Ca — fio as function of 
the saddle's order fraction Us/SN for MLJ (open symbols) 
and BMLJ (full symbols). The dashed line, with a slope of 
~ 11 is the best fit to the MLJ data, b) The mean square 
distance between minima (IS) reached by quenching the 
inherent saddle of order Us is plotted as function of ns/3N 
for MLJ. The dashed line is the best fit to the data 
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